Nonadiabatic corrections to the wave function and energy 
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Nonadiabatic corrections in molecules composed of a few atoms are considered. It is demon- 
strated that a systematic perturbative expansion around the adiabatic solution is possible, with the 
expansion parameter being the electron- nucleus mass ratio to the power 3/4. Closed form formulae 
for the leading corrections to the wave function and to the energy are derived. Their applicability 
is demonstrated by a comparison of numerical results for the hydrogen molecule with the former 
nonadiabatic calculations and the experimental values. Good agreement with the recent experiment 
is achieved for the ground state dissociation energy of both H2 and D2 . 
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I. INTRODUCTION 



One of the most fundamental approximations of the 
theory of molecular structure, the so caUed clamped 
nuclei approximation, assumes separation of the elec- 
tronic and nuclear motion and computing the electronic 
wave function and energy for a fixed position of nuclei 
P, [3, H, 01 • In the more accurate adiabatic approxima- 
tion, one includes the diagonal matrix element of the nu- 
clear kinetic energy, see Eq. ([8|). The reminder beyond 
this term is described as the nonadiabatic correction. In 
their calculations Kolos and Wolniewicz Q pointed out 
significance of nonadiabatic corrections in obtaining ac- 
curate results for H2 . They claimed that within the Born- 
Oppenheimer (BO) approximation or adiabatic variants 
thereof, it is not possible to estimate quantitatively the 
accuracy of the method, and suggested that one should 
drop from the very beginning the idea of separation of 
the electronic and nuclear motions, and start instead 
with the exact Schrodinger equation for all the parti- 
cles involved. Since that work, calculations including the 
nonadiabatic as well as the relativistic effects in the hy- 
drogen molecule have been significantly im pro ved, among 
others, by Wolniewicz and Kolos [ESBBISL Adamow- 
icz [lllIlllilH, Bishop [11, Anderson [11 and their 
coworkers. However, no systematic theory of nonadia- 
batic corrections to dissociation energies or rovibrational 
splittings has been pursued. For simple molecules such as 
H2, direct nonadiabatic calculations with ex plic itly cor- 
related Gaussians are at present possible [lll.ll2l. n3l . [l^ , 
but no reliably estimate of their numerical uncertainties 
has been made. In fact, results from Refs. [l2| and [l3| 
are not in good agreement with each other, if all pre- 
sented digits are significant. However^ the recent results 
of the variational calculation in Ref. [ij are in remark- 
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able agreement with the original work of Wolniewicz [3] , 
but our results obtained here for total nonrelativistic en- 
ergies lie slightly below their predictions, outside the es- 
timated uncertainties. 



At present, the lack of nonadiabatic corrections in sim- 
ple molecules limits the knowledge of physical properties, 
such as the dissociation energy, vibrational, and rota- 
tional energy spacings. In certain instances the nona- 
diabatic effects can not be considered small at all. For 
example, a recent study of the spontaneous ortho-para 
transition in H2 [l3| reveals that the contribution of the 
nonadiabatic effects to the transition probability reaches 
as much as 40%. Another example concerns the helium 
dimer. The nonadiabatic corrections are expected [11 
to significantly contribute to the dissociation energy of 
the highest vibrational level D14 of the spin polarized 
helium molecule in electronic state, for which precise 
measurement has recently been performed [l9| . Indeed, 
authors of ^23] observed that using the atomic mass in- 
stead of the nuclear mass in Eq. ^ j,] changes the dis- 
sociation energy of the D14 level by as much as 5%. This 
mass change is in accordance with the result obtained 
here, Eq. (|33p . and indicates an importance of the nona- 
diabatic effects. 



In this work, we construct a systematic perturbative 
expansion for nonadiabatic corrections to the wave func- 
tion and energy in molecules composed of a few atoms, 
with the expansion parameter being the electron-nucleus 
mass ratio to the power 3/4 [{nic/M)^^^]. We derive 
closed form formulae for the leading nonadiabatic cor- 
rections and present accurate numerical results for the 
H2 and D2 molecules. This method, we suppose, can be 
applied to more complex molecules, where fully nonadi- 
abatic calculations are infeasible, which is the main ad- 
vantage of the presented here perturbative approach. 
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II. METHODOLOGY 
A. Theoretical framework 

Let us first consider a diatomic molecule, and assume 
the reference frame be related to the mass center of the 
two nuclei. The total wave function is the solution of 
the stationary Schrodinger equation 



with the Hamiltonian 



H — Hc\ + ffn : 



(1) 



(2) 



split into the electronic and nuclear parts. In the elec- 
tronic Hamiltonian iJoi 



E 



2 Top 



V 



(3) 



the nuclei have fixed positions Ra and Rb, while the 
nuclear Hamiltonian is 



Hn = - 



V2 



V2 



2 Ma 2 Mb 



2^„ 2 {Ma 



Mf. 



(4) 
(5) 



where R = Ra — Rb and /in is the nuclear reduced mass. 

In the adiabatic approximation the total wave function 
of the molecule 



^,{r,R)^M^RX{R) 



(6) 



is represented as a product of the electronic wave function 
(/>ci and the nuclear wave function x- The electronic wave 
function obeys the clamped nuclei electronic Schrodinger 
equation 



[i/el - Scl(i?)] |0cl) = 0, 



(7) 



while the nuclear wave function obeys the Schrodinger 
equation in the effective potential generated by electrons 



(8) 



To account for corrections to the adiabatic wave function 
0a and the energy E^. resulting from the coupling be- 
tween the motion of nuclei and electrons, we formulate a 
systematic perturbative expansion around the adiabatic 
solution, and derive explicit expression for the leading 
nonadiabatic corrections to the wave function 5(/)na and 
to the energy 5E^g_. 



B. Derivation of tlie nonadiabatic corrections 



The total wave function 

= 0a + 54>ria. = 0cl X + ' 



(9) 



is the sum of the adiabatic solution and a nonadiabatic 
correction. Due to the normalization of (j) and 0aj the 
nonadiabatic correction is orthogonal to 0a in the leading 
order of the electron-nucleus mass ratio 



W„a|0a) =0. 

Let us consider the electronic matrix element 



^nal^^cl/cl 



(10) 



(11) 



which is a function of the internuclear distance R. The 
nonadiabatic correction (50na will be decomposed into two 
parts 



= 0cl 5X + ^'011 



(12) 



It follows from Eqs. (fTO|) and (fTT|) that these two parts 
obey the following orthogonality conditions 



(l5'0na|0cl)cl = 0, 

{5x\x) = 0. 



(13) 
(14) 



Let us now revert to the stationary Schrodinger equation 
in (HI) 



[H^l + i7„ - Sa - 5E^^] |0cl X + '50„a) - , 

and rewrite it to the following form 

[E,,{R)-H4 |^0„a) 

= [E,x{R) + i/„ -E,- 5E^,] |0ei X + ' 
The solution can be formally written as 

1 



(15) 



(16) 



|<50na) = |0ol'5x) 



[EcI + Hn — Es. — SEn 



(17) 



|0elX + ^0na), 



where the symbol prime in the denominator denotes ex- 
clusion of the reference state 0ei from the Hamiltonian 
inversion. Since (50na is expected to be a small correction 
to 0a, it can be neglected on the right hand side in the 
leading order. Then Ed — Eg^ — 5 Ens, does not contribute 
in the above matrix element and this equation takes the 
form 



|<50na) = |0C1<^X) + 



1 



{E, 



-ff„|0eix). (18) 



Let us point out, that the difference E^i — E^, in Eq. pT)) 
is not necessarily small for higher rovibrational levels, so 
the neglect of i50na on the right hand side might be ques- 
tionable. However, we will show later, that this neglect 
of (50na is in fact justified by considering the next or- 
der nonadiabatic correction. The presence of Ec\ instead 
of ii'a in the denominator in Eq. (|18p is an important 
difference with the former calculations of nonadiabatic 
corrections by Wolniewicz [7]. 
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The perturbed function, Sx, is as yet unknown — we 
obtain it from Eq. (jlSp by taking a matrix element with 
(pel, namely 

{c^,l\E,i + Hn-E,-SEn4 0cl (x + '5x)+'5'0na),i = (19) 

), which can be 



The nuclear function x obeys Eq. ([8] 
rewritten in the form 

{4'oi\EcI + i?n - £^a| 4>c\ x) ^ 



0. 



(20) 



The difference of Eqs. (fT9|) and ([20ll . neglecting the small 
higher order term SE^a Sx is 

[i7n + (i/„),i + i?cl-i?a] |<5x) (21) 
= SEnaX ~ (0cl|-ffn|<^'0na)^i- 

If one multiplies the above equation from the left by x 
and integrates over i?, then obtains 



(5-Ena = (0cix|i?n|(5'0na) 



1 



(E^\ — Hoi) 



001 X - (22) 



This is the final formula for the leading nonadiabatic cor- 
rection to energy. At first sight, it resembles the Van 
Vleck sum-over-states formula 0, [Hj, but the essential 
difference is in that, no summation over the nuclear states 
is present here, which makes Eq. ([22]) of practical use. 
The perturbed function 6x can now be inferred from 
Eq. (HH). The solution to this differential equation can 
be formally written as 



I^X) = 



[Ea — Eel — Hn — (-ffn)pi] 



-{(t)cl\Hn\S'(pna) 



(23) 



This completes the treatment of leading nonadiabatic 
corrections in diatomic molecules. The extension of 
the perturbative nonadiabatic expansion to polyatomic 
molecules is straightforward. Let the total Hamiltonian 
be decomposed into the electronic and nuclear parts as 
in Eq. ([2|), with the condition that -ffoi does not contain 
derivatives with respect to nuclear coordinates. Then the 
nonadiabatic correction to energy again takes the form 
of Eq. ((221) • This general formula is also valid in the case 
when the i?-independent term in in Eq. ([5]) is shifted 
to Hoi, which can simplify the numerical calculation of 
nonadiabatic energies. 



C. Higher order corrections 

The perturbative theory presented here can be ex- 
tended also to higher orders in the electron-nucleus mass 
ratio in a manner similar to the standard time indepen- 
dent Rayleigh-Schrodinger perturbative expansion. In- 



deed, the next order nonadiabatic correction to the en- 
ergy is 



S^^^Ena = (</'clX 



1 



{Eel — Hel) 



7 (^^n + Ed — Ea) 



{Eel — Hcl) 



0olX 



0clX 



H 



" {Eel — Hel 



(beiSx), (24) 



where Sx is given in Eq. (|23p . and its detailed derivation 
is postponed to the next paper. Let us note, that H^ + 
Eei^Ea in the first term of the right hand side of Eq. ((24|) 
is indeed small, because x satisfies differential equation 
([5]). Therefore, the neglect of S4)na on the right hand side 
of Eq. ([T7)l in the derivation of the leading nonadiabatic 
correction is fully justified. 

The order of magnitude of the higher order corrections 
can be estimated by noting that in atomic units the ki- 
netic energy of the nuclei is of order of the square of the 
nuclear displacement SR from the equilibrium value. It 
follows that SR ~ (mc/Ain)^^"', and the nuclear kinetic 
energy is nominally of order y^mT/A^n- However, it can 
be shown that each H^ in Eqs. ((22|l and ((24|l contributes 
(TOo/Mn)^^^, therefore the nonadiabatic correction to the 
energy is of order ©[(rTioZ/in)"^^^] , and the total energy 
can be written as 



E^Ea + SEn 



l + O 



3/4 



(25) 



This estimate may break down however, when Born- 
Oppenheimer potential curves for different electronic 
states approach each other at certain distance, then the 
denominator in Eq. (|22p can be arbitrarily small. For 
this situation a kind of multireference perturbation the- 
ory should be developed and this is not a subject of the 
present work. 



D. Working equations 



The general formula ([22|) can be readily rearranged to 
a more practical form. For this purpose, we separate out 
electronic matrix elements from the nuclear ones in SE^a 
and assume that the states x are purely vibrational 



SE„a - 



i?-^ dR 



1 



Mn \ 

v^xv|x/™ 



{Eel — Hel 
1 



Hn I 



\ 



{Eel - Hex)' 
1 



V|j0el ) (26) 
cl 



V^<?!>ci 



[Eel - Hel)' 

i?2 dR [x^R) U{R) - 2 x(i?) x'(^) V{R) 



[X'{R)?W{R)]. (27) 
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Eq. ((27)) can be rewritten as an integral of the sum of two 
terms which have simple physical interpretation 



+ [X\R)YW{R) 



U{R) + -V[R)+V{R) 



(28) 



The first term, the coefhcient of x^, is a nonadiabatic 
correction to the adiabatic energy curve, while the second 
term, the coefficient of x'^, is an i?-dependent correction 
to the nuclear reduced mass /in, namely 



1 



1 



2/in(i?) 2/i„ 



W{R). 



(29) 



An obvious advantage of the representations ([27|) and 
([28]) is that the electronic matrix elements comprising 
the pseudopotentials U, V, and W need to be evaluated 
for each R only once for all rovibrational levels. 

The electronic matrix elements in Eq. (j26p involve mul- 
tiple differentiation of the electronic wave function with 
respect to the internuclear distance i?, which is difficult 
to calculate directly. Therefore, we rewrite these terms 
to a more convenient form, where differentiation is taken 
of the Hamiltonian, namely 



1 



RVcl 



{Eel — Hel 



(30) 

- I V^ (ilel) (/-cl + 2 Vfl(i/el - £^cl) 
-Vi^(iIel)0cl|, (31) 



1 

{Ed — Hci) 
1 



where, in the last equation, a term proportional to (/)ci 
has been omitted as it does not contribute to the matrix 
element in Eq. ([26]) . 



E. Asymptotics 

The adiabatic energy Es, in Eq. ([8]) and the nonadia- 
batic correction SEm, in Eq. do not vanish at large 
internuclear distances. It is because they partially in- 
clude the atomic reduced mass and the so called mass 
polarization term. For example, for the large atomic sep- 
aration in the hydrogen molecule, (i?n)ci in Eq. ([5]) and 
U{R) in Eq. (|26p are equal to mc/mp and — (mc/'rip)^, 
accordingly, which corresponds to first terms in the ex- 
pansion of the reduced mass in the electron-nucleus mass 
ratio 



1-J^ = 

rric 1 -|- mc/mp m^ 



mc/mp _ mo _ / ^"\^ ^ f _ 
^mp) \mp) 

(32) 

The asymptotics of W{R) equal to —mg/mp corresponds 
to the change in Eq. (|29p of the reduced nuclear mass /!„ 



to the reduced atomic mass /ia, namely 
1 1 1 1 



1 - 



2 /iA 2/in(oo) mp + 771c " 

(33) 

Therefore, if one uses, instead of the nuclear reduced 
mass the atomic reduced mass fj,A in the radial equa- 
tion for X, then the large R asymptotics should be sub- 
tracted from W{R). 



III. CALCULATIONS 

In order to test validity of the presented perturba- 
tion theory, we consider in more detail the hydrogen 
and deuterium molecules, for which the nonadiabatic 
corrections were calculated perturbatively by Wolniewicz 
Til, and nonperturbatively by Adamowicz and coworkers 
13, [13, [ll El. The nuclear masses used in our calcu- 
lations were mp = 1836.15267247 rrie for the proton and 
TOd = 3670.4829654 m^ for the deuteron, and the conver- 
sion factor was 1 hartree = 219474.6313705 cm~^ 

The electronic potential for the nuclear Schrodinger 
equation ([8]) consists of two contributions: the BO en- 
ergy curve Eel and the adiabatic correction curve (i/n)ci- 
In our calculations we used the most accurate cur- 
rently available BO potential of Ho computed with sub- 
nanohartree accuracy by Cencek [23| using 1200-term ex- 
jxmentially correlated Gaussian (EGG) wave functions 

The adiabatic corrections for H2 and D2 were com- 
puted as expectation values of the Hamiltonian Hn, 
Eq. ©, 



{4>cl\Hn\(j)cl)c 



-^(0cl|(Vi+V2)'|(/.el)el. (34) 



While the evaluation of the last term is straightforward, 
the first term on the right hand side is cumbersome since 
it involves differentiation of the electronic wave function 
with respect to the internuclear distance. Several meth- 
ods of evaluation, either analytically 0, |2^ [2^ or nu- 
merically i28i, |29,] of the adiabatic correction (i?n)ci 
have been reported in literature. In this work, we note 
that the expectation value (0ci|V^|0ei)ci, by virtue of 

Vfl((/)oi|</>ci)oi = 0, is equivalent to -(Vfl^^oil Vfl,(/)ei)ei, and 
with the help of Eq. (|30p is expressed as 



l|Vfl(i7el)' 



1 



[{Eel^Hcl)'] 
xVfl(iJel) |0cl), (35) 



which can be conveniently calculated in a similar way as 
the nonadiabatic correction. 

The numerical calculations were performed at a few 
tens of internuclear distances, R. For each R, three dif- 
ferent electronic functions expanded in an EGG basis set 
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were involved in the evaluation of the components of the 
vector Vfl</)ei of Eq. (PH)). The two-electron ECG basis 
functions were of the following form 



Mn,r2) = (i + Pi2)(i + J)Sfc 

2 

X exp 



(36) 



where the matrices and vectors Sk contain nonlinear 
parameters, 5 per basis function, to be optimized. The 
antisymmetry projector {I + P12) ensures singlet symme- 
try, the spatial projector (1 -I- i) — the gerade symmetry, 
and Sfc prefactor enforces S state when equal to 1, and 
n state when equal to yi — the perpendicular Cartesian 
component of the electron coordinate. 

The first basis, composed of 600 ECG functions ([36|) . 
was employed to expand the X"'^E+ electronic ground 
state wave function c/fd. The nonlinear parameters were 
optimized variationally with respect to the BO energy 
with the target accuracy of the order of nanohartree. 
Two other basis sets were employed to represent the 
two components of Vfl(/)ci: a 1200-term basis set of 
symmetry for the component parallel to the molecular 
axis and a 1200-term basis set of ^Ilg symmetry for the 
perpendicular component. The nonlinear parameters of 
these two basis sets were also optimized. This optimiza- 
tion was carried out with respect to the following func- 
tional 



1 



- Vr{H,{) |(/.el) (37) 



with the fixed external wave function (/)ci and in the case 
of the S basis set, with the subtraction of the internal 
wave function 0oi from the Hamiltonian inversion. This 



is achieved as follows. The first 600 terms of this basis 
were taken from (pi^i wave function and their nonlinear 
parameters were kept fixed during the optimization, only 
the remaining 600 terms of E symmetry were optimized. 
This ensures that the internal wave function (pd is well 
represented at every step of optimization. Then the sub- 
traction is achieved by orthogonalization of V;? (i?ci) l^oi) 
with respect to the internal \4>c\)- In the final calculations 
the 1200-term E basis obtained according to the above 
description was used also for the expansion of the exter- 
nal ground state function <f>c\- 

The optimized functions were employed to evaluate 
the adiabatic correction ([34|) yielding an excellent agree- 
ment with the most accurate results obtained by Cencek 
and Kutzelnigg |29j by using the numerical Born-Handy 
method — all six digits reported in [2^ were confirmed by 
our analytic calculations. The diagonal correction func- 
tion obtained this way was combined with the BO energy 
curve to form the electronic potential for the movement of 
the nuclei. The nuclear Schrodinger equation was solved 
numerically for the adiabatic energies E^, and the nuclear 
wave functions x(^)- 

The electronic matrix elements iY, V, W entering 
Eq. ((27)) were evaluated within the same two 1200-term 
ECG basis sets described above yielding smooth func- 
tions of R. Because for the highest vibrational levels the 
nuclear wave functions are spread out and contributions 
from larger internuclear distances are non-negligible, the 
functions U{R), V{R), and yV(i?) were represented by 
their asymptotic forms: U{R) = uq + uq/R^ + us/R^, 
V{R) = ve/R^ + vs/R^, and W{R) ^ wq + we/R^ + 
ws/R^, subject to uq = Wo = —{nic/M)'^ restriction (in 
atomic units). The remaining, free parameters u;, w;, 
and Wi were determined by fitting the above functions to 
calculated points in the range of (5.5 : 10) bohrs. 
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IV. RESULTS AND DISCUSSION 

Results of the computations are collected in Tables |T] 
andini where the nonadiabatic corrections SEm^ obtained 
for all purely vibrational levels of H2 and D2 are listed 
and compared with the existing data. We observe small 
deviation from results of Wolniewicz Q for both H2 and 
D2. Apart from the highest vibrational levels, the devi- 
ation is about —0.013 cm~^ for H2 and —0.003 cm~^ for 
D2. It can be attributed to the slightly different pertur- 
bation approach of Wolniewicz in Ref. [3]. Since, the fully 
nonadiabatic calculation by Stanke et al. in Ref. is 
in remarkable agreement with results of Ref. for H2, 
the most probable explanation of this deviation is the ne- 
glect of the probably significant second term in the higher 
correction of Eq. ((24)) . We note, that with the varia- 
tional approach the obtained results are upper bounds 
to the nonadiabatic energies and the estimate of uncer- 



tainty is not necessarily simple, as shows the comparison 
with the former calculations by Bubin and Adamowicz 
in [1^. While we cannot judge here which result is more 
accurate, we observe, after inclusion of relativistic and 
QED corrections [sO] , a very good agreement of our final 
values with the precise experimental result for the ground 
state dissociation energy in both, H2, see Table [llD and 
D2 molecule in Table IIVI 



A. Simplified approach 

In the case of the lowest rovibrational states, the nona- 
diabatic correction of Eq. ()26p can be significantly sim- 
plified. Since the derivatives of <j)c\ are here much smaller 
than that of x, the contribution from terms proportional 
to xx' and x^ can be neglected. Moreover, the nuclear 
wave function x is strongly localized around i?0 7 the av- 
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TABLE I; Comparison of the nonadiabatic corrections to the adiabatic energy levels of H2, computed using Eq. (|22|) . with 
results of Wolniewicz [3] as well as of Stanke et al. (in cm~^). Relative uncertainty of our results due to the neglect of 
higher order corrections is 0.36%. In the second column, our reference adiabatic energy _Ea is given for all the vibrational levels. 



V 




faul 


L'-'^iia 


f7l 


Diff. 


fl4l 


Diff. 





_1 




—0.511 


—0.4988 


—0.012 


—0.499 


—0.012 




_1 


yjtjij ZiOU 


— 1.347 


— 1.3350 


—0.012 


— 1.336 


—0.012 


2 


— 1 


127168 403 


—2.104 


—2.0913 


—0.013 


—2.092 


—0.012 


3 




1 1 327 S41 


—2.785 


—2.7728 


—0.013 


—2.773 


—0.012 


A 
ft 


1 

— _L 


OQzL P19^ 7^17 


— o.oyo 




u.u 10 


— 0.000 


n ni 9 


5 


-1 


079 751 576 


-3.934 


-3.9208 


-0.013 


-3.921 


-0.013 


6 


-1 


066 017 251 


-4.398 


-4.3847 


-0.013 


-4.385 


-0.013 


7 


-1 


053 339 044 


-4.779 


-4.7654 


-0.013 


-4.765 


-0.013 


8 


-1 


041 750 041 


-5.059 


-5.0459 


-0.013 


-5.046 


-0.013 


9 


-1 


031 301 694 


-5.210 


-5.1980 


-0.012 


-5.197 


-0.013 


10 


-1 


022 068 803 


-5.187 


-5.1772 


-0.010 


-5.175 


-0.012 


11 


-1 


014156 653 


-4.926 


-4.9176 


-0.008 


-4.915 


-0.010 


12 


-1 


007 711399 


-4.333 


-4.3270 


-0.006 


-4.324 


-0.009 


13 


-1 


002 935 436 


-3.287 


-3.2844 


-0.002 


-3.281 


-0.006 


14 


-1 


000 108 439 


-1.646 


-1.6482 


0.002 


-1.645 


-0.001 



TABLE II: Comparison of the nonadiabatic corrections to the 
adiabatic energy levels of D2, computed using Eq. (|22p . with 
results of Wolniewicz 0] (in cm~^). Relative uncertainty of 
our results due to the neglect of higher order corrections is 
0.21%. In the second column, our reference adiabatic energy 
is given for all the vibrational levels. 



lying vibrational states is 



V 


Sa[au] 


D2 

(5-Bna 


\^\ 


Diff. 





-1.1671680230 


-0.176 


-0.1725 


-0.003 


1 


-1.1535267225 


-0.480 


-0.4769 


-0.003 


2 


-1.1404282296 


-0.764 


-0.7605 


-0.003 


3 


-1.1278630079 


-1.028 


-1.0246 


-0.003 


4 


-1.1158235203 


-1.273 


-1.2698 


-0.003 


5 


-1.1043043280 


-1.500 


-1.4970 


-0.003 


6 


-1.0933022295 


-1.710 


-1.7062 


-0.004 


7 


-1.0828164445 


-1.901 


-1.8977 


-0.004 


8 


-1.0728488533 


-2.075 


-2.0709 


-0.004 


9 


-1.0634043039 


-2.228 


-2.2246 


-0.004 


10 


-1.0544910050 


-2.361 


-2.3568 


-0.004 


11 


-1.0461210239 


-2.468 


-2.4646 


-0.004 


12 


-1.0383109190 


-2.547 


-2.5434 


-0.004 


13 


-1.0310825465 


-2.590 


-2.5869 


-0.004 


14 


-1.0244640945 


-2.590 


-2.5868 


-0.003 


15 


-1.0184914179 


-2.534 


-2.5316 


-0.003 


16 


-1.0132097815 


-2.410 


-2.4071 


-0.002 


17 


-1.0086761569 


-2.197 


-2.1951 


-0.002 


18 


-1.0049622689 


-1.875 


-1.8742 


-0.001 


19 


-1.0021585658 


-1.420 


-1.4196 


0.000 


20 


-1.0003781627 


-0.805 


-0.8055 


0.000 


21 


-0.9997348811 


-0.079 


-0.0860 


0.007 



erage distance between nuclei, so the electronic matrix 
element can be taken just at i?o- The resulting simpli- 
fied form of nonadiabatic correction to energy for low 



h^n \ 



[{E,, - Fel)'] 



Tl- Vfl(i?ol) 

3 n ■ Vfl(i?ol) 



yd ; , (38) 



where 
n-Vfl(i7ci) 



n f + (39) 

^3 ' ^3 ^3 ^^3 ] jj)2 ' \ ' 



lA 'IB ' 2A ' 2B 



and n = R/R. The corresponding result, obtained for 
the ground vibrational state of H2, —0.453 cm"-'^ differs 
from the accurate one (with the asymptotics subtracted 
out) —0.4458 cm~^ only by 1.5%. Moreover, this approx- 
imation can be applied also to heavier molecules, where 
its accuracy should be even higher. Indeed, the simplified 
calculations for the ground state of D2 yield —0.161 cm^^ 
which differs from the accurate result —0.1594 cm~^ by 
as few as 1.0%. 



V. SUMMARY 

We have formulated the perturbative approach to the 
calculation of nonadiabatic effects. The leading order 
nonadiabatic correction to the energy, given by Eq. (|22p . 
is conveniently rewritten in terms of three state indepen- 
dent pseudopotentials introduced in Eq. (|27p . and cal- 
culated in Sec. IIIII for all vibrational levels of H2 and 
D2. Our numerical results are slightly below, and thus in 
small disagreement with calculations by Wolniewicz Q 
and by Stanke et al. [l^ . 

Several important advantages of the approach pre- 
sented here are worth emphasizing. The first one is the 
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possibility to systematically derive higher order nonadi- 
abatic corrections and to control its accuracy. Secondly, 
the result for the leading correction of Eq. ([22|) can conve- 
niently be applied to more complex molecules. Thirdly, 
for the lowest vibrational states one can apply the ap- 
proximate Eq. ([38]) . which is accurate to about 99%, even 
for very light molecules such as H2 and D2. Finally, the 
nonadiabatic corrections to the wave function Eq. (fT5)) 
significantly contribute to radiative transition rates, such 
as ortho-para transition in H2 [l7j . or make the transi- 
tion possible, as in the case of rovibrational electric dipole 
transitions in HD. 
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TABLE IV: Dissociation energy of D2 u = state (in cm~^). 
Total uncertainty is due to higher order nonadiabatic cor- 
rections (0.0004 cm~^) and approximate a* QED correction 
(0.0008 cm-i). 

Correction Subtotal 



BO 

+ adiab. corr. + asymp. 
+ nonad. corr., Eq. 
4- asymp. (|IIEp 
+ o? relat. corr. [30] 
o? QED corr. [30] 
+ a* QED corr. [30] 
uncertainty 

Experiment [31] 
difference 



0.1757 
-0.0163 
-0.5178 
-0.2359 
-0.0017 



36746.162 

36748.935 

36749.111 

36749.094" 

36748.577 

36748.341'' 

36748.339 



±0.0012 36748.339(1) 

36748.343(10) 
0.004(11) 



"Bubin and Adamowicz 'X1\: 36749.085 

''Wolniewicz [7| (includes approximate QED corr.): 36748.364 

in Poznah Supercomputing and Networking Center. 



TABLE III: Dissociation energy of II2 = state (in cm~^). 
Total uncertainty is due to higher order nonadiabatic cor- 
rections (0.0019 cm^^) and approximate a** QED correction 
(0.0008 cm-i). 

Correction Subtotal 



BO 



-I- ad. corr. -(- asymp. E|) 
-I- nonad. corr., Eq. 
-I- asymp. E|) 
-f o? relat. corr. [3^ 
-f o? QED corr. [30] 



QED corr. [30] 



0.5109 
-0.0651 
-0.5178 
-0.2359 
-0.0017 



36112.593 

36118.364 

36118.875 

36118.810" 

36118.292 

36118.056'' 

36118.054 



uncertamty 



±0.0027 36118.054(3) 



Experiment [31] 
difference 



36118.062(10) 
0.008(13) 



"Bubin and Adamowicz 'l2^: 36118.798 

''Wolniewicz 0| (includes approximate QED corr.): 36118.069 
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